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We show that the dynamics of interacting fermions can be exactly replaced by a quantum jump 
theory in the many-body density matrix space. In this theory, jumps occur between densities formed 
of pairs of Slater determinants, D a t = |$a) where each state evolves according to the Stochastic 
Schrodinger Equation (SSE) given in ref. Q|. A stochastic Liouville-von Neumann equation is 
derived as well as the associated Bogolyubov-Born-Green-Kirwood-Yvon (BBGKY) hierarchy. Due 
to the specific form of the many-body density along the path, the presented theory is equivalent 
to a stochastic theory in one-body density matrix space, in which each density matrix evolves 
according to its own mean field augmented by a one-body noise. Guided by the exact reformulation, 
a stochastic mean field dynamics valid in the weak coupling approximation is proposed. This theory 
leads to an approximate treatment of two-body effects similar to the extended Time-Dependent 
Hartree-Fock (Extended TDHF) scheme. In this stochastic mean field dynamics, statistical mixing 
can be directly considered and jumps occur on a coarse-grained time scale. Accordingly, numerical 
effort is expected to be significantly reduced for applications. 
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I. INTRODUCTION 



The purpose of this paper is to discuss the possibility to substitute the description of the evolution of quantum 
interacting fermions by a stochastic mean-field dynamics of one-body density matrices. In view of present computa- 
tional capabilities, stochastic methods appear as a promising tool to address exactly or approximately the problem 
of correlated mesoscopic quantum systems such as nuclei, atomic clusters or Bose-Einstein condensates. Mean field 
theories, i.e. Hartree-Fock theories, are rarely able to describe the large variety of phenomena occurring in quantum 
systems. It is generally necessary to extend mean field theory by including the effect of two-body correlations 0]. 
During thepast decades, several approximate stochastic theories have been proposed to describe strongly interacting 
systems 0, 0, H, El H> El EH [HI- These approaches have in common that the noise is due to the residual part of 
the interactions acting on top of the mean field. However they generally differ on the strategy used to incorporate 
noise. In some cases, the residual interaction is treated using statistical assumptions 0,0, while in other cases the 
interaction induces fluctuations in the wave- packets either by random phase-shifts or b^quantum jumps according 
to the Fermi-Golden rule 0, ^| . The influence of correlations is sometimes treated using the notion of stochastic 
trajectories in the one-body density matrix space0, 0- This latter is, among the different theories, the only one 
that has been applied to large amplitude collective motions in the semi-classical limit [l^. Recently, its quantal 
version has been used to describe small amplitude collective vibrations in nuclei |l3[ . However the application of a 
stochastic approach to the quantum many-body dynamics remains an open problem both from a numerical as well as 
a conceptual point of view jlj, [l5j . 

In this work a different strategy is used to obtain a stochastic formulation of the many- body problem. During the 
last ten years, many efforts have been made using functional integral techniques 0,0,03 to address the problem of 
nucleons in strong two-body interactions. These theories provide an exact stochastic formulation of quantum problems 
and lead to the so-called quantum Monte-Carlo methods Recent applications to nuclear physics have shown that 
stochastic methods can be applied successfully to describe the structure of nuclei [2(j. These methods can also be 
applied to the description of dynamical properties 0] . However the self-consistent mean field does not generally play 
a special role. Indeed the stochastic part is driven either by the kinetic energy part of the Hamiltonian or by a fixed 
one-body potential in the case of shell model Monte-Carlo calculations [2jj. Recently a new formulation has 
been proposed that combines the advantages of the Monte-Carlo methods and of the mean field theories. Application 
of functional integral theories are of great interest since they pave the way to a full implementation of the nuclear 
static and dynamical many-body problem using mean field theories in a well defined theoretical framework. However 
the direct application of exact stochastic dynamics to realistic situations remains numerically impossible and proper 
approximations should be developed. 
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In the first part of the article, the functional integral method and the associated Stochastic Schroedinger Equation 
(SSE) developed in Q for many-body and one-body wave functions are presented. The theory is formulated in the 
more general framework of exact stochastic dynamics in the many-body and/or in the one-body density matrix space. 
The link between the different formulations is underlined. In a second part, guided by the exact stochastic theory, an 
extended mean field theory |14| taking into account two-body effects in the weak coupling regime is given in terms of 
a new stochastic one-body evolution. 



II. INTRODUCTION AND DISCUSSION OF STOCHASTIC METHODS 



Functional integral methods have been used for a long time to provide a useful reformulation of complex quantum 
system s flfil IT^ | (for a review see 0])- This method has been applied with success to describe static properties of 
nuclei 20]. However, it was seldom used for dynamical problems. Recently, an alternative formulation of the path 
integral representation has been obtained in which the mean field theory plays a specific role. We consider a general 
many-body system described by the wave function |<I>) which evolves according to the Hamiltonian 1 : 

H = Y^ T ii a t a 3 + V ljk ia+a+aia k , (1) 

ij ijkl 

where the first term corresponds to the kinetic part of the Hamiltonian while the second part is the two-body 
interaction. We use the convention of [22| concerning the labelling of one and two-body operators. We denote 
Vijki — (ij \vi2 1 kl) where v\2 is the antisymmetrized two-body interaction. 



A. Action of a quadratic Hamiltonian on a Slater determinant 

In ref. [2(j, the general strategy to obtain ground state properties of a many-body system using Monte-Carlo 
methods is described. The new aspect developed in ref. is the introduction a self-consistent mean field before the 
application of functional integral. In that case, only the residual part of the interaction which is not taken into account 
in mean field Hamiltonian is treated stochastically. In this section, we summarize how a general two-body Hamiltonian 
applied to a Slater determinant can be separated into a mean field part and a residual two-body contribution. Details 
are given in ref. [J]. 

We consider a Slater determinant |$) defined as 1$) = n Q a+ |0) , where the single particle states \a) may not be 
orthogonal. Starting from the Hamiltonian we have 

H\&) =(Hi+H res )\<f>), (2) 

with 

Him =\E +Y, (ai\h MF (pi)\a 1 )a+ 1 a &1 |$), (3) 

\ QiaT / 

where we denote by the particle states (i.e. the unoccupied states) and where pi —^2 (6ti\ is the one-body 
density associated with |$). The states |di) are defined by (ct\ \ a.2) = 5 aia . 2 . In this expression, Hmf (pi) is the 
mean field Hamiltonian 

hMF(pi) = Tt +v{pi) . (4) 

In this equation, v (pi) = TT2 (V12P2) is the mean field potential where Tr2 (.) denotes the partial trace on the second 
particle. In equation J3J, we denote 

E a = Tr (pih MF (p x ) - l -piv( Pl )) . (5) 



Note that three-body (or higher) interactions are not considered here. 
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In the single particle basis defined above, it could be shown that 

(ai<5 2 |«i 2 | aia 2 ) = - ^ tiWa (ai \O s \ ai) (a 2 |O s | a 2 ) , (6) 

s,aict2c*ic*2 

where O s is a one-body operator 2 . Note that, the latter transformation of the two-body matrix elements is a particular 
case of the more general transformations given in ref. |2(J . When the single particle basis is not the particle-hole state 
of the Slater determinant, additional terms should be accounted for. Using this transformation, the residual part of 
the Hamiltonian is 

Hres \®) = \ huJs ( Sl ai ) a 2> a^a^as!^ |$) . (7) 

In next section, this expression is the starting point to derive the stochastic Schroedinger equations using functional 
integral techniques. 



B. Functional integrals and stochastic many-body dynamics 

Functional integrals methods applied to quantum fermionic systems in interaction |l6l ITt) lead to general stochastic 
formulations of the quantum many-body problem. They however also lead to specific difficulties. For instance, the 
semi-classical limit of the functional integral does not give naturally the Hartree-Fock theory, but only to the Hartree 
theory. The interesting idea proposed in [ij is to use the functional integral already accounting for the fact that 
the Hamiltonian is applied to a Slater determinant. In this case, only the residual (2 particle-2 hole) part of the 
Hamiltonian is interpreted as a source of noise. This procedure is summarized now. 

We consider the evolution of the system during a small time-step At. Denoting by |A$) the associated evolution, 
we have: 

|$ + A$) = U(At)\&) =cxpf^-H) |$), (8) 



where U (At) is the propagator associated to H. Using the Hubbard-Stratonovitch |23|,|24!] functional integral on the 
residual part only, the exact propagator transforms into an integral equation Q] : 

U (At) |$) = / dlfG (if) exp (f H x + AB (if)) |*>. (9) 

Hi is given by equation Q while AB (if) is a one-body operator written as 

AB(o>)|$) = X * (a\O s \a) AW s ata & \<5>) , (10) 

where 

\ s = ^/uT s (l + isgn(uj s )) /2 (11) 

and AW S — vAt a s with a s the component of the vector T? . In equation G (if) = H s g(cr s ) represents the product 
of normalized Gaussian probabilities of width 1 for the a s variables. As in other functional integral formulations, we 
recover that the original propagator associated to the exact evolution can be replaced by an ensemble of propagators 
that depend on if. Equivalently, in the limit of infinitesimal time step (At — ► dt), this equation can be interpreted 
as a stochastic Schroedinger equation for the initial state. Using the standard notation for stochastic processes in 
Hilbert space [25|, we have: 

|$) + |d$> = exp (^H x + dB (t)) |$) . (12) 



n 



2 Following ref. |22| . we will sometimes make use of the identity U12 = — ^ HuigOlO 2 which is a compact notation for matrix elements 
and is only valid in the particle-hole basis. Here, we use the same notations as in ref. Il-ll where "1" and "2" refers to the particles on 
which the operator is acting. 
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Here \d$>) has to be interpreted as a stochastic wave function. Since eq. I© is exact, it shows that the exact dynamics 
of a Slater determinant can be replaced by an average over stochastic evolution operators. In this expression, dB (t) 
is a stochastic operator which depends on the stochastic variable dW, according to equation (|lf)|) 3 . In order to obtain 
this equation, the Ito rules for stochastic calculus have been used [2f| with 

dW Sl dW S2 = 6 slS2 dt. (13) 

Using the latter properties in combination with the expression of dB (t) , we obtain an equivalent of the fluctuation- 
dissipation theorem that gives the link between the stochastic operator and the residual part of the Hamiltonian: 

\dB (t) dB (t) |*) = +%H res (t) |$) . (14) 
2 in 

Expression (|12fl is of particular interest. Indeed, according to the Thouless theorem [26l |2?|. the application of 
an operator of the form (|12|) to a Slater determinant gives another Slater determinant. Therefore, the evolution of 
correlated systems of fermions can be replaced by stochastic evolutions of an ensemble of Slater determinants. Since 
each evolution can be solved with numerical techniques used in mean field theories, SSE offer a chance to solve exactly 
the dynamics of strongly interacting fermionic systems. This property has already been noted in several pioneering 
works [Ifil IT?! ITsj . A very similar conclusion has been reached for the description of interacting bosons using Monte- 
Carlo wave function techniques [2l| • In this case and more generally in the context of the stochastic description of 
open quantum systems, jumps between wave-packets are generally described using differential stochastic dynamics in 
Hilbert space 28, 29, 30]. Then, the evolution of \d<S>) is directly considered. 

The equivalent differential equation associated to the jump process described here can be obtained by developing 
the exponential in eq. (|12f> in powers of dt. Using Ito rules, we obtain : 

|*) + |d$) = ( 1 + %H X + \dB (t) dB (t) + dB {t)) |$) . (15) 
\ in 2 J 

Using equations and l|14fl . we finally obtain a stochastic Schroedinger equation for the many-body wave function: 

|d*> = (jt H + dB(t)J |$). (16) 

In the following, this equation is referred to as the many-body SSE. Equation (|16[) is strictly equivalent to i|12|) and thus 
preserves the Slater determinant nature of the states along the stochastic trajectory. This might appear surprising due 
to the appearance of the complete Hamiltonian H in eq. ill6|) . This is a specific aspect of the stochastic many-body 
theory using Ito stochastic calculus. Indeed, although H (which contains the complete two-body interaction) drives 
the initial state out from the Slater determinant space, the stochastic part of the equation of evolution compensates 
this effect exactly. The exponential form l|12fl or the differential form 116fl describe the same stochastic process. 
However, differential equations are generally easier to manipulate p5 l29l l3fjf . 



C. Equivalent quantum jump for single particle states 



Up to now, we have introduced notions associated with the stochastic mechanics of many-body wave functions. 
This formulation is of great interest for applications since the stochastic evolution of the many-body wave function can 
be replaced by the stochastic evolution of its single particle components. For completeness the equivalent differential 
equation of single particle wave function is given below. It has been shown in ref. |l( that equation (|12fl leads to the 
single particle equation of motion 

|da) = \ a ) + K (1 - Pi) O s \a) dW s , (17) 

s 

where h (pi) is a one-body operator given by 

h(px) = h M F{p\)-\piv{p\). (18) 



3 Note that in the limit At » dt, dW s plays directly the role of the Gaussian normalized stochastic variable, and the introduction of a s 

is not required. 
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Equation l|17|) will be referred to as the one-body SSE. We would like to stress again that eq. I|16[) and the set of single 
particle evolutions (eq. (|17|l ) are strictly equivalent. 

In this section, we have summarized the equivalence between quantum jump approaches in many-body and one- 
body spaces of wave-packets in order to describe interacting fermions. An equivalent formulation in terms of density 
matrices is highly desirable to compare the exact treatment with other stochastic methods. 



III. DENSITY MATRIX FORMULATION 



In the previous section, we have considered the stochastic formulation of the many-body problem using stochastic 
Schroedinger equation. In this approach, all trajectories start from a Slater determinant and follow a stochastic path 
in the Slater determinant space. Stochastic theories can also be applied if the system is initially correlated. In this 
case, it is helpful to generalize the theory by introducing density matrices. It has been shown in ref. that the 
many-body density matrix D (t) associated with the system at all times can be properly described by the average 
over an ensemble of pairs of non-orthogonal Slater determinants state vectors, 

D(t)=WJWl (19) 

each of them evolving according to eq. (|16() . Here, the average over the initial ensemble has been introduced. In 
that case, the notion of a quantum jump between the wave functions is replaced by a quantum jump in the space of 
Slater determinants pairs. In the following, the properties of Slater-determinant dyadics are recalled and a stochastic 
BBGKY hierarchy iHIHE! is derived. 



A. Slater determinant dyadics: notations 



Let us consider a many-body density formed of two distinct Slater determinants 

Dab = |$.) ($&| , (20) 

in which each Slater determinant is an antisymmctrized product of not necessarily orthogonal single particle states 

|$ a ) =n afl +|o) 

|$ 6 ) =n,a+\o) ■ W 

Note that Z? a {, is neither hermitian nor normalized. However, for convenience we will still call it a density matrix. 
Starting from the many-body density matrix, one can obtain the generalized k— body density matrix (denoted by 
Pi...k ) by taking successive partial traces. Using the same notation as in j^], we have: 

P ( " b i=Tr k+1 ..., A (D ab ), (22) 

where A is the size of the system. One can obtain the expression of density matrices in terms of single particle states 
of the two Slater determinants by introducing the overlap matrix elements between single particle states, denoted by 
/. The matrix elements of / are defined by f/3 iaj = For instance, the one-body density matrix is: 

P ^ b) = det (/) £ la,) (ft, | = det (/) u^ b) . (23) 

ai/3j 

More generally the k— body density matrix is the antisymmetrized product of single particle densities |34| 

P i ^lk = tet(f)A(u[ ab) ...u < £V), (24) 

where A(.) corresponds to the antisymmctrization operator. Introducing the two-body correlation operator defined 
by 

C^=p^-A{u^p^), (25) 
we have = for any state defined by equation (|20|l . 
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B. Stochastic evolution of many-body density matrices 

The BBGKY hierarchy 0, I32I l33| has been widely used as a starting point in order to obtain approximations 
[2^ on the evolution of complex systems. Therefore, an equivalent hierarchy associated to the exact stochastic mean 
field deduced from functional integrals is highly desirable to specify the possible links with other theories. In this 
section, starting from the stochastic Schroedinger equation for the many-body wave function, we give the associated 
stochastic formulation of the BBGKY hierarchy. In the stochastic many-body dynamics, we consider the quantum 
jump between two different density matrices D a b and D' ab . Starting from D ab given by eq. (12011 . there are transitions 
towards another density matrix given by D' ab = |$ a + rf$ a ) ($& + The rules for transitions are directly obtained 

from the rules for the jumps in the wave functions space: 



</<!>„} = #/i |<i>„> +<ib„ <]>„ 

with 



(26) 



(d* 6 | = -§ {<S> b \H + {$ b \dB+. 
dB a =Y,s X sT l &aai(a\Os\a}a & dW Sa 

dB+ = T, s K »J (P Pb\ P) apdW Sb . ^ 

The notation dW Sa and dW Sb are introduced in order to emphasize that stochastic variables associated respectively 
to |$ a ) and are statistically independent, i.e. 

dW s JW Sb = 0. (28) 

This complete eq. Ijl3|) verified both by dW Sa and dW Sb . With these rules, the evolution of the many-body density 
matrix along the stochastic path is given by 

dD ab = ^ [H, D ab ] + dB a D ab + D ab dB+. (29) 

This equation is a stochastic version of the Liouville-von Neumann equation for the density matrix. The evolution of 
the k— body density matrix can be directly derived from expression (|29|) and one obtains: 

d Pi b i = ^Tr k+1 ,..., A ([H, D ab ]) + dW k ab . (30) 

The additional term corresponds to the stochastic part acting on the fc-body density matrix evolution: 

dW k ab = Tr k+1 _ A (dB a D ab + D ab dB+) . (31) 

The first part of eq. I|30|) is nothing but the standard expression of the k th equation of the BBGKY hierarchy whose 
explicit form can be found in review articles p2l l35l l36| . The equation of motion for the fc— body density matrix in 
the framework of the stochastic many-body theory proposed in ref. [l| corresponds to the standard BBGKY term 
augmented by a one-body stochastic noise. 



C. Evolution of the one-body density matrix 



Starting from (|3L)fl . an explicit form of the one-body density evolution can be found. Since for any D ab , we have 

ab 
12 



Ci2^ = 0' tne nrs t term in eq. l{3*0"j) reduces to 



Tr 2 ..., A ([H,D ab \) =Tr 2 



H,p[ a 2 b) 



(ab)\ (ab) 



Pi 



(32) 



The stochastic part reads: 



Tr 2 ,..., A (dB a D ab ) = E S A S £ C 



10* 



'Tr 2 ,...,A (ata & D ab ) dW Si 



(33) 
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Let us introduce a complete single particle basis. For any state \i) and \j) of the basis, we have : 

(i \Tr 2 ,....A (ata & D ab ) \ j) = Tr (a+ aia%a & D ab ) 

= J2ki ( k I "> (" I Tr {a+a.a+aiDab) . 



(34) 



Using the fermionic commutation rules on creation/annihilation operators together with the definition of the one- and 
two-body density matrices, we obtain: 



(i \Tr 2 



.,A [a 



ta&Dab) | j) = J2h (k\a)(a\l)(l 



I) (i 



(ab) 
Pl2 



(ab) 
Pi 



(35) 



Using the fact that p\i = A (uxp2), we finally obtain: 



Tr 2 ,..., A (dB a D ab ) = E s A s ( 1 - u[ ab) ) 0.tf°>dW t 



■£ S A S TY (u[ ab) (l- Pa )O s )p[ ab) dW Si 



(36) 



where the one-body density p a associated with |$ a ) has been introduced. The same treatment can be performed for 
the second part of the stochastic term and the evolution of the one-body density matrix finally reads: 



dp\ 



(ab) 



dt 

if 7 ! 



T-MF 



b) 



with 



db 



(ab) _ 



tab) 



oJ* b) dw l 



(u[^),p[« 



db 



(ab) 



(ab) 



(l- Pa )O s )p[ ab) dW t 



E s Kp[ ab) O s (l - <4° 6) ) dW Sb + E s KTr (O, (1 - pb) u[ ab) ) p[ ab) dW Sl 



(37) 



(38) 



It is interesting to note that although the single particle states entering in p^^ do not evolve according to mean 

field theory but according to h {^>i b ^ given by l|18l) . the deterministic part associated with the evolution of the 

one-body density reduces to the standard mean field propagation. Eq. I|37|) points out the central role played by the 
mean field Hamiltonian in the stochastic many-body theory. In particular, it shows that any evolution of a correlated 
physical system submitted to a two-body interaction can be replaced by a set of mean field evolutions augmented by 
a one-body noise. Finally, it is worth noticing that expression l|37|l can alternatively be obtained by differentiating 

directly p{ ab) = det (/) J2 a ^ \at) (&|. 



D. k-body density evolution from one-body density 

The stochastic evolution transforms a pair of Slater determinants into another pair of Slater determinants. Thus, 
all the information on a single stochastic trajectory is contained in the stochastic evolution of the one-body density 
evolution in eq. (|37l) . Indeed, the evolution of the fc— body density matrix can be directly obtained from the relation 
(I24|l which is valid all along the stochastic path. Using the Ito rules, we have 

dp^ k = d(det(f))A(u[ ab K..u[^) 

+dct( f )j: l A(u^ b) ...du^...ut b) ) 

+d(dct(f))J2 i A(u^...du < f b \..u{ ab ^ 

+ det (/) A (E^ u[ ab) . . . du[ ab) . . . duf } . . . u[ ab) ) . 

It can be checked that the terms which are linear in dt correspond to the deterministic part of eq. (|30|l . The latter 
expression is also useful in order to have an explicit form of the stochastic noise to all order in k. In expression 139fl . 
d(det (/)) is deduced from equations l|2t)|) . We have 

ddet(f) = (d^> b | $„) + ($ 6 1 d$ a ) (40) 
= (® b \dB a + dB+\$ a ), (41) 
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Many-Body 
Stochastic Schrodinger equation 



Stochastic evolution 
of many-body density 



dD n 



m 



[S, D ab \ + dB a D ab 4- D ab dB+ 



|tfe) = %h. u>\) w + £ -v, (i - pi) o s .«> m, 

in g 

One-Body — 
Stochastic Schrodinger equation 



Stochastic evolution 
of one-body density 



P \ ab n+db^ 



FIG. 1: Summary of the four different ways of considering the exact reformulation of the quantum many-body problem using 
stochastic mechanics. The single arrow indicates that the density matrix formulation can be derived from the stochastic 
Schroedinger equations. The double arrows show that for both wave function or density matrix formulations, there is a strict 
equivalence between the many-body and the one-body stochastic equation of motion. 



which gives 



ddct (/) = X *Tr (u ( f b) (1 - Pa ) O a ) dW Sa + KTr (o s (1 - p b ) u[ ab) ) dW s , 

S S 

In addition, the equation on du[ is deduced from (|37f) . Altogether, we obtain 



(42) 



1 - u^ ab) ) O 



dW s . 



+ Es ^Tr (u[ ab) (1 - Pa ) OA dt^ Sa pi' 
+ EsKTr (OAl~p b )u { « b) )dW SbP [ a 



(Ob) 



(43) 



Here, we introduced the notation 0\ to denote that the one-body operator O s is applied to the particle "i". The 
possibility to derive the evolution of p\...k for all k from the evolution of p\ is an illustration of an attractive aspect 
of this theory. Indeed, since we are considering pairs of Slater determinants, all the information on the dynamics 
is contained in their one-body densities. This proves that the exact evolution of the density matrix of a correlated 
system through a two-body Hamiltonian can always be replaced by the average evolution of uncorrelated states each 
of them evolving in the one-body space according to its own mean field augmented by a one-body stochastic noise. 



E. Summary 

Functional integral methods are attractive since they provide a rather transparent and systematic way of trans- 
forming the exact dynamics of a correlated system into a stochastic mean field dynamics. In this work, we have 
discussed the link between the different one-body and many-body SSE's on one side and the stochastic one-body and 
many-body density evolution on the other side. The equivalence and the relationship between the various ways of 
considering stochastic mechanics are displayed in fig. ^ 

The exact stochastic formulation of the dynamics of complex systems provides a well defined framework to introduce 
stochastic theories. However, the stochastic dynamics as it is proposed is still rather cumbersome as far as numerical 
applications are concerned. Indeed, due to the increasing number of trajectories with the number of degrees of freedom, 
exact stochastic many-body theories have only been applied to dynamics of rather schematic models !]. With present 
computational facilities, there is no chance to apply the exact theory to realistic mesoscopic systems and approximate 
formulations are necessary. The stochastic theory provides however a natural way to replace the dynamics of an 
interacting system by one-body dynamical evolutions. In the following, we will transform the stochastic equation to 
account approximately for the correlation and reduce the numerical effort. 
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IV. APPROXIMATE STOCHASTIC MANY-BODY DYNAMICS 



A number of approximations of the many-body problem can be found in the literature. Among them, the mean 
field theory is certainly the most widely used. Correlations beyond the mean field are often required to have a realistic 
description of dissipative aspects in mesoscopic systems. A general strategy to obtain extensions of the mean field 
dynamics consists in performing successive truncations of the BBGKY hierarchv |l4ll35ll37| . The first order truncation 
of the hierarchy leads for instance to the standard mean field theory. An extension of the mean field can be obtained 
by considering the first and second equations of the hierarchy. This has led to different levels of approximations to 
the nuclear many-body problem as for instance the so-called extended Time-Dependent Hartree-Fock[3|||3!j^(J (f° r 
a recent review see ref. [lj])- In the following, we will show that the stochastic evolution described previously can be 
adapted to a stochastic one-body theory for correlated systems equivalent to the extended TDHF. 



A. Extended mean field dynamics 

Theories beyond mean field [HEl 

are valid when the dynamical effect of the residual interaction is weak. In the 
weak coupling regime, correlations can be treated perturbatively on top of the mean field. These theories are valid 
under the assumption that different time scales associated respectively to two-body collisions and to the mean field 
propagation exist. Consider t co u the time scale for an in- medium two-body collision, and Tf ree the time between two 
collisions. In the weak coupling approximation, one can assume that there exists a time interval At verifying the 
condition 

Tcoll < At < Tfree- (44) 

An estimate and a discussion of these time scales can be found in ref. 0, 0] . The physical picture to interpret 
the separation of time scales is that each single particle state evolves according to the average mean field and rarely 
"encounters" a two-body collision. From the many-body problem point of view, the role of the residual part of the 
interaction is to account for two-body collisions. 

Besides time scales, extended TDHF remains a one-body theory. Indeed, it is assumed that part of the two- 
body correlations can be neglected and that the two-body density matrix can be instantaneously approximated by 
an antisymmetrized product of one-body density matrices (puit) — A [pi(t)p2 (*)])• This is of special interest for 
practical applications since only one-body degrees of freedom are followed in time. 



B. Approximate stochastic dynamics 

In this section, we propose a formulation of extended one-body dynamics in terms of quantum jumps in the space 
of one-body density using the same hypothesis as in extended TDHF. We start from a system described at time to 
by its one-body density given by 

Pi {to) = i a ) n « H • ( 45 ) 

a 

The system is assumed to be initially uncorrelated so that p\2 (to) = -A (pi (to) pi (to))- Let us now consider an 
ensemble of one-body density matrices, noted p± with initial conditions p\^ (t§) — p\ (to). The time interval At is 
divided into N time steps (At = N As) and at each time step, p^ evolves according to its mean field augmented by 
a stochastic term 

Ap^ = ^[h MF (p^),pM}+AK(pM). (46) 

However, contrary to the strategy of the previous section, and following the hypothesis of extended mean field theory, 
jumps are supposed to occur only once in the time interval At. For a jump occurring at a time r = fcAs, the stochastic 
term is written as 



AK (p[ n) (t)) = 8 t , T {E s A s (l - p[ n) (t)) O s pi n) (t) AW S 

+ J2s'K'p[ n) (t)o s (i-p[ n) (t))&w s ,}, 



(47) 
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where A s and O s are denned in previous section, while AW S and AW S > are two independent Gaussian stochastic 
variables that follow Ito stochastic rules, with 



AW S1 AW S2 = S SlS2 As 



AW,' AW,, 



■As. 



(48) 



We consider the ensemble of trajectories with a quantum jump occurring at a specific time r. pJ 2 W denotes the 
two-body density obtained by averaging A p™ p% over these trajectories. 

Before the time r, all trajectories follow the same path corresponding to the mean field propagation with the initial 
condition p\ (to). We note respectively p™^ an d U m f the associated one-body density and propagator. We have 



P? f (t') = U mf (f, t ) Pl (t ) U+ f (f, to) , 



with 



ds 



(")>') 
Pi 



U m f (t',t ) = Texp (^+~^ /i ™/(p™ / ( S )) 

Using these definitions, the evolution between r and r + As of the product (p^' 

A (piV) = (ApW) P ^ + A n) (Aa4 b#) ) + (Api" } ) (Ap^" 
Using expression (|46[) and Ito rules, we obtain 

ApfW) = f [ft m/ (pr / w)+ftm/(^ / w),pr / wpr / w 

+AX (p( n) ) AX (p| n,) ) + AX (p<">) + AK (p { f> 



(49) 



(50) 



(51) 



(52) 



We have used the fact that, for all considered trajectories, no collision occurs before time r leading to Pi(t) — 

P\ ^( T ) — pT* ( t )- The last two terms of equation (|52"|l do not contribute to the average evolution. We thus see that, 
in addition to the mean field, an extra deterministic term will appear in the average evolution l|52(l . Using equations 
\, we have 



AK (p[ n) ) AK (p 2 n ">) = AsJ2 s (\ s y(l-p?f(T)) (l-pT f (r))0l0 2 sP T f (r) P T f (r) 



A S E s (K) 2 P? f (r) PT (r) 0\Ol ( 1 - p™> ( T ) ) ( 1 - ft* (t 



"if 



(53) 



Using finally the fact that = iu> s /2, relation (jSJ) and introducing the antisymmetrization operators, we obtain the 
average evolution 



As 
if i 



hmf ( P T f (r)) + h mf (p^(r)) ,A(pT f (T)p^(T)) 



AA(p^ P ^)(t) 

In this equation, Fi 2 reads 

Fi2(r) = (l - p[ nf (r)) (l - P 7 f (r)) v 12 p^ (r) ft* (r) 
-P? 1 (r) P 2 f (r) v 12 (l - P T f (r)) (l - (r)) 



As 



Fu (r) . 



(54) 



(55) 



As discussed in [l5|, the effect of a single collision is expected to be weak during the time interval At and we can 
assume that for all trajectories, the mean field propagation coincides with U m f after the jump. Therefore, the average 
density at the final time tf = to + At is given by: 



P12 



(t f ) = A (p? f (t f )p 2 nf (tf)) + ^„ 2 / (*/,r) F 12 (r) U%+ (t f , r) , 



(56) 



where = U m j ® U m j. The complete average density p\ 2 (tf) is obtained by summing different possible times r 
for collisions: 



(*/) = A (p? f (tf)pf(t f j) + ^ f S dsU™ (t s , s) F 12 (s) U%? (t f ,s) , 



Pl2{ 



(57) 
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where the limit As — > ds has been taken. This two-body density matrix corresponds to the standard mean field 
propagation augmented by the incoherent contribution of nucleon-nucleon collisions entering generally in extended 
mean field theories 

As mentioned previously, an interesting aspect of extended TDHF is that it contains only one-body degrees of 
freedom. This can only be achieved by projecting correlation effects in the single particle space. In the stochastic 
dynamics presented here, this is equivalent to assume that the final two-body density can be approximated by 
•A(pi(tf)p2(tf)) where pi(t f ) is given by 

Pl (t f )=Tr 2 (m(tf)). (58) 

The density obtained in this way differs from the density propagated by the mean field alone and contains the effect 
of incoherent nucleon-nucleon collisions. The procedure can then be iterated using the new density as a starting point 
for future stochastic propagation. 

In this section, we have presented a method to include approximately two-body effects by means of a stochastic one- 
body theory. As in the exact formulation presented in the previous section, the stochastic theory can be equivalently 
formulated as a stochastic Schroedinger equation. It is important to note that the numerical effort required for the 
approximate dynamics is expected to be much less than for the exact one at least for two reasons. The first one comes 
from the fact that quantum jumps occur on a "coarse-grained" time-scale. The second reason lies in the possibility of 
directly propagating densities formed by a statistical mixing (eq. I@SJ) without invoking pairs of Slater determinants. 
As a counterpart, we would like to mention that the approximate stochastic formulation has the same limitations at 
extended TDHF and can only be applied to problems for which the residual correlations are weak. 

V. CONCLUSION 

The main result of our work is the proof that the exact dynamics of a correlated system evolving through a two- 
body Hamiltonian can be replaced by a set of stochastic evolutions of one-body density matrices where each density 
evolves according to its own mean field augmented by a one-body noise. Guided by the exact stochastic formulation, 
an approximate stochastic mean field theory valid in the weak coupling limit is proposed. In this theory, jumps occur 
on a coarse-grained time scale. 

The alternative stochastic formulation presented here does avoid some of the ambiguities present in other stochastic 
theories. A first remarkable aspect comes from the justification of the noise source. Indeed, since the starting point 
of our work is an exact formulation of the many-body problem, the noise has an unambiguous mathematical and 
physical interpretation. 

In addition, from a practical point of view, it has clearly some advantages. In all applications to quantum problems 
of extended mean field theory, it has been shown that the memory effect is important (see discussion in |l4l Il5jl. 
This memory effect corresponds to the non-local action in time of the past history collisions on the future dynamics. 
Although the noise is Markovian, it accounts also for this non-Markovian effect through the instantaneous average 
over trajectories. In addition, as noted in ref. 0, in order to apply stochastic theories proposed in ref. @, H(| to 
large amplitude motions, one should be able to guess what will be the important states in the future evolution. This is 
in particular necessary to reduce the number of trajectories. For instance, it has been guessed in ref. |ld| that jumps 
can be optimized due to the 2 particle-2 hole (2p-2h) nature of the residual interaction. In the theory developed here, 
the system is driven naturally towards the important states. Indeed, as can be seen from eq. ()17fl . these states are 
self-consistently defined without ambiguity, and the 2p-2h character of the residual interaction directly shows up in 
the stochastic part of the propagator. 

The exact treatment of the many-body problem with stochastic theories is still not possible for realistic large 
amplitude dynamics due to the required numerical effort. However, an alternative formulation of the stochastic 
theory has been proposed in the second part of this article which should make the numerical applications easier. This 
stochastic theory provides a suitable framework for the description of interacting systems in the weak coupling regime. 
In particular, it keeps the advantages discussed above and it is expected to significantly reduce the numerical efforts 
for practical applications. Such a theory could a priori be applied to nuclear systems where quantum and dissipative 
effects are important such as for instance giant resonances, fusion reactions or the thermalization in nuclear reactions. 

Finally we would like to mention that an additional difficulty may be encountered due to the possible progressive 
entanglement of the initial state. Indeed, starting from an initial simple state, the states propagated with stochastic 
Schroedinger equation will progressively become more complicated and fragmented over phase space. If such an 
entanglement occurs, the method proposed here might be very difficult or even impossible to use. 
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